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Abstract 

We develop the cosmological perturbations formalism in models 
with a single non-local scalar field originating from the string field 
theory description of the rolling tachyon dynamics. We construct the 
equation for the energy density perturbations of the non-local scalar 
field in the presence of the arbitrary potential and formulate the lo- 
cal system of equations for perturbations in the linearized model when 
both simple and double roots of the characteristic equation are present. 
We carry out the general analysis related to the curvature and entropy 
perturbations and consider the most specific example of perturbations 
when important quantities in the model become complex. 



1 Introduction 



Current cosmological observational data [IJ [2] strongly support that the 
present Universe exhibits an accelerated expansion providing thereby an ev- 
idence for a dominating dark energy (DE) component [3]. Recent results 
of WMAP [2] together with the data on la supernovae and galaxy clus- 
ters measurements, give the following bounds for the DE state parameter 
wde = — 1-02^q j 6 . The present cosmological observations do not exclude a 
possibility that the DE with w < — 1 exists, as well as an evolving DE state 
parameter w. Moreover, the recent analysis of the observation data indi- 
cates that the varying in time dark energy with the state parameter wde, 
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which crosses the cosmological constant barrier, gives a better fit than a 
cosmological constant (for details see reviews [1] and references therein). 

Construction of a stable model with w < — 1 is a challenge leading to 
the consideration of originally stable theories admitting the NEC violation 
in some limits. Recently a new class of cosmological models obeying this 
property which is based on the string field theory (SFT) [5] and the p-adic 
string theory [6] has been investigated a lot [7]-[26]. It is known that both 
the SFT and the p-adic string theory are UV-complete ones. Thus one can 
expect that resulting (effective) models should be free of pathologies. 

Models originating from the SFT are distinguished by presence of specific 
non-local operators. The higher derivative terms in principle may produce 
the well known Ostrogradski instability |27j (see also [281 [8])0. However 
the Ostrogradski result is related to higher than two but a finite number 
of derivatives. In the case of infinitely many derivatives it is possible that 
instabilities do not appear [19]. 

The SFT inspired cosmological models [7] are considered as models for 
dark energy (DE). The way of solving the Friedmann equations with a 
quadratic potential, by reducing them to the Friedmann equations with 
many free massive local scalar fields, has been proposed in [10[ [T2] (see 
also [25]). The obtained local fields satisfy the second order linear differen- 
tial equations. In the representation of many scalar fields some of them are 
normal and some of them are phantom (ghost) ones [121 125j . Cosmological 
models coming out from the SFT or the p-adic string theory are considered 
in application to inflation [l5]-[22] to explain in particular appearance of 
non-gaussianities. For a more general discussion on the string cosmology 
and coming out of string theory theoretical explanations of the cosmologi- 
cal observational data the reader is referred to |31| . Other models obeying 
non- locality and their cosmological consequences are considered in [32\ [33] . 

As a simplest model originating from SFT one can consider a theory with 
one scalar field whose kinetic operator is non-local. Equations for cosmolog- 
ical perturbations in such kind of model where the scalar field Lagrangian 
is quadratic covariantly coupled with Einstein gravity were derived in |23j . 
In the present paper we develop and improve that formalism accounting an 
arbitrary potential of the scalar field as well as the presence of double roots 
of the characteristic equation in the linearized model. We also carry out the 
general analysis of curvature and entropy perturbations and consider the 
most specific example of perturbations when characteristic quantities of the 
model become complex. 

The paper is organized as follows. In Section 2 we describe the non-local 
non-linear SFT inspired cosmological model. In Section 3 we sketch the con- 
struction of background solutions in the linearized model and perturbation 

Additional phantom solutions, obtained by the Ostrogradski method in some models 
can be interpreted as non-physical ones, because the terms with higher-order derivatives 
are regarded as corrections essential only at small energies below the physical cutoff [291130] . 
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theory for models with non-local scalar fielco In Section 4 we consider the 
perturbations in the case of two complex conjugate roots. In Section 5 we 
consider the case when the analytic function ^(c/D), which defines the ki- 
netic operator in the action, has double roots. In Section 6 we summarize 
the obtained results and propose directions for further investigations. 



2 Model setup 
2.1 Actions 

We work in (1 + 3) dimensions, the coordinates are denoted by Greek indices 
[A, v, . . . running from to 3. Spatial indexes are a,b, . . . and they run from 
1 to 3. The four-dimensional action motivated by the string field theory is 
as follows [351 dOl E] : 

s =\*~l=i (isb + 7ft {¥ r " (a ' n)T - V ^ T) ) - A ») • (1) 

Here Gat is the Newtonian constant: 8ttGn = 1/Mp, where Mp is the Planck 
mass, a' is the string length squared, g a is the open string coupling constant. 
We use the signature (—,+,+,+), g^ v is the metric tensor, R is the scalar 
curvature, □ = D^d^ = —^d^^—gg^dy and being a covariant deriva- 
tive, T is a scalar field primarily associated with the open string tachyon. 
The function J-Q^a'O) may be not a polynomial manifestly producing thereby 
the non-locality. Fields are dimensionless while [g ] = length. The poten- 
tial Vi n t(T), which is an open string tachyon self- interaction, does not have 
a quadratic term. It is convenient to introduce dimensionless coordinates 
= x^/yfa! , the dimensionless gravitational constant Gn = Gn/ol , and 
the dimensionless coupling constant g Q = g j\fal ' . In the following formulae 
we always use dimensionless coordinates and parameters omitting bars over 
them. 

Function Tq is assumed to be an analytic function of its argument, such 
that one can represent it by the convergent series expansion with real coef- 
ficients: 

oo 

Jo = ^/nD n and f n G R. (2) 
Equations of motion are 



n=0 



_ 8itG n 

9o 

T (D)T = VUT), (4) 



2 For applications of other multi-field cosmological models and related technical aspects 
see for instance 1341. 
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where is the Einstein tensor, the energy-momentum (stress) tensor is 
as follows 



n-l 



71=1 1 = 

- g»v (g^d p U l Td a U n - l - l T + u l TU n ~ l T^ 

- g»» (&Ao - Qrj- (D)r - v M (T))\ . 

It is easy to check that the Bianchi identity is satisfied on-shell and for 
Fq = /iD + fo the usual energy-momentum tensor for the massive scalar 
field is reproduced. 

Let us emphasize that the potential of the field T is V = — ^-T 2 +Vi n t{T) . 
Let To be an extremum of the potential V. One can linearize the theory in its 
neighborhood using T = Tq + r. To second order in r one gets the following 
action 

where F = F G - V(T )" and A = A + Equations © and Q are 

9o 

valid for the latter action after the replacement J-q — > F and Ao — > A at 
Vi n t(T) = 0. Note that all Taylor series coefficients f n , except fo, are the 
same for Fo and F. Equation @ now is 

F(D)t = 0. (6) 

Non-local cosmological models of type (|5|) with 

F sit (a) = -fa + i- ce ~ 2n , 

were previously considered in [TTJ [121 H7pl . Actions (HJ and ([5]) are of the 
main concern of the present paper. 



2.2 Background solutions construction in the linearized mo- 
del 

While solution construction in the full non-linear model ([T]) is not yet known 
the classical solutions to equations coming out the linearized action © were 
studied and analyzed in ITJJI HU H3 HZl I2D1 1231 Thus, we just briefly 
notice the key points useful for purposes of the present paper. 

3 In [17] for example it has been shown that solving the non-local equations using the 
localization technique is fully equivalent to reformulating the problem using the diffusion- 
like partial differential equations. One can fix the initial data for the partial differential 
equation, using the initial data of the special local fields. This specifies initial data for 
a non-linear model, and these initial data can be (numerically) evolved into the full non- 
linear regime using the diffusion-like partial differential equation. 
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The main idea of finding solutions to the equations of motion is to start 
with equation ([6]) and to solve it, assuming the function r is a sum of eigen- 
functions of the d'Alembert operator: 

r = T i> where Drj = JiTi and ^{Ji) = for any i = 1, . . . , N. (7) 

i 

Hereafter we use N (which can be infinite as well) denoting the number of 
roots and omit it in writing explicit summation limits over i. Without loss 
of generality we assume that for any i\ and i% ^ i\ condition J- n ^ Ji 2 is 
satisfied. In this way of solving all the information is extracted from the roots 
of the characteristic equation J 7 (J) = 0. We can consider the solution r as a 
general solution if all roots of J- are simple. The analysis is more complicated 
in the case of double roots [25]. We consider this case separately in Section 
5. 

In an arbitrary metric the energy-momentum tensor in (|3|) evaluated on 
such a solution is [23] 

Tpv = ?(Ji) (fynOvn - 9 -f {g pa d p T t d a Ti + J t T 2 )) - glg^A. (8) 

i 

The last formula is exactly the energy-momentum tensor of many free mas- 
sive scalar fields. If J-{ J) has simple real roots, then positive and negative 
values of J-'(Ji) alternate, so we can obtain phantom fields. Using formula 
(jSJ) we obtain the Ostrogradski representation [271 EE] for action (0): 

53 =/ ^ (idb " A " g W (ra ' TATi + J - Tf >) ■ 

one can see that S3 is a local action if the number of roots N is finite. 

2.3 Application to Priedmann— Robertson— Walker Universe 

We stress that all the above formulae are valid for an arbitrary metric and the 
general solution. From now on, however, the only metric we will be interested 
in is the spatially flat Friedmann-Robertson-Walker (FRW) metric with the 
interval given by 

ds 2 = -dt 2 + a 2 (t) [dx\ + dx\ + dxf) (9) 

where a(t) is the scale factor, t is the cosmic time. 

Background solutions for r are taken to be space-homogeneous. The 
energy-momentum tensor in ([3|) in this metric can be written in the form of 
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a perfect fluid Ty = diag(— Q,p,p,p), where 

-. oo n—1 

q = -^2f n J2( d t^ l Td t a n - 1 - l T + n l Ta n - l T^j - 

n=l 1=0 

- l -TT Q {U)T + V mt {T) + g 2 o A , 

2 (10) 

OO 71—1 

p = - f n (d t U l Td t U n - 1 - l T - U l TU n ~ l T\ + 

n=l 1=0 

+ ^TFo{U)T -V int (T) - g 2 o A . 

Obviously we can rewrite equations ([3]) in the canonical form: 

3H 2 = 8ttGq, H= -4ttG{q + p), (11) 

where G = G^jg 2 is a dimensionless analog of the Newtonian constant, 
H = a I a and a dot denotes a derivative with respect to the cosmic time t. 
The consequence of (jlip is the conservation equation: 

g + 3H(g + p) = 0. (12) 

which corresponds to the non-local field equation @. 

Some progress is possible in the linearized model when the metric is fixed 
to be of the FRW type. The individual equations in ([7]) in this metric read 

n + 3Hf t + J in = (13) 

The full system of equations of motion has the fixed points at n = and 
3H 2 = 3Hl = 87tGatA, which is real at A > 0. Equations $13]) together 
with Friedmann equations describe the late time evolution of the model 
with Lagrangian ([T]). This model possesses an approximate solution with 
all scalar fields tending to the minimum of the potential (i.e. t% — > 0) and 
the Hubble parameter going to the constant. Such solution was constructed 
numerically and was proven to be a solution in p3]. Also the asymptotic 
form of this solution was derived in [101. 

It is instructive to investigate the Lyapunov stability of the fixed point. 
Using formulae from [4.1j, we come to conclusion the fixed point is asymp- 
totically stable at 

H > 0, sfte(Ji) < 0. (14) 

If at the fixed point Hq < or 3?e( Jj) > for some i, then this fixed point is 
unstable. At ffie(Ji) = for some i or Hq = one can not use the Lyapunov 
theorem to analyse the stability of the fixed point. Note that the conditions 
(|14p are sufficient for stability in not only the Friedmann-Robertson-Walker 
metric but also the Bianchi I metric |41| . In this paper we shall extend this 
analysis and consider the stability of the fixed point with respect to arbitrary 
perturbations. 
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To compute an approximate solution to (|13p one starts with a constant 
H = Hq and then computes the correction to H using Friedmann equations. 
Then the procedure can be iterated to compute higher corrections. It was 
proven in [10] such iteration does converge. 

Solution to (fT3"j) with constant H = Hq is obviously 

n = r i+ e Q «+* + Ti-e**-* (15) 



where a«± = ^M 3 - (— 1 ± A /l — ). Considering r, we see that only one 



2 

term in the solution converges when t — > oo in general (if both converge we 
select the slowest one). Let's assume it is the first one proportional to r + 
constant. Then in order to pick the (slowest) convergent solution we put 
r_=O0 

The first correction to the constant Hubble parameter in case only de- 
caying modes in n are present is 

H = H + h = H + h J2 T i- ( 16 ) 

Constant ho is not independent and is related with n+. We note that h 
is of order rf . The last expression is a good approximation for H in the 
asymptotic regime when h <C Hq. Further one can find the scale factor to 
be 2 

a = a exp(W+yE^:)- ( 17 ) 

3 Cosmological perturbations with single non-local 
scalar field 

3.1 General analysis 

Now we turn to the main problem of the present paper: derivation of cosmo- 
logical perturbation equations in models with a non-local scalar field. We are 
focused on the scalar perturbations, because both vector and tensor pertur- 
bations exhibit no instabilities [37J . Scalar metric perturbations are given by 
four arbitrary scalar functions [361 [57] . Changing the coordinate system one 
can both produce fictitious perturbations and remove real ones. Natural way 
to distinct real and fictitious perturbations is introducing gauge-invariant 
variables, which are free of these complications and are equal to zero for a 
system without perturbations. There exist two independent gauge- invariant 
variables (the Bardeen potentials), which fully determine the scalar pertur- 
bations of the metric tensor [3J3 [37J [381 HOI [23] . To construct the perturba- 
tion equations one can use the longitudinal (conformal-Newtonian) gauge, 



4 Hereafter we adopt the rule \fz* = y/z* meaning that the phase of the complex number 
runs in the interval [— 7T, n) and for z = re !CT the square root is |y / r|e !CT/ ' 2 . 
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in which the interval Q with scalar perturbations has the following form 
(in terms of the Bardeen potentials): 

ds 2 = a( V ) 2 (-(1 + 2$)d V 2 + <^(1 - 2^)dx a dx b ^ (18) 

where r] is the conformal time related to the cosmic one as a(r])dr] = dt. The 
the Bardeen potentials and are as usually Fourier transformed with 
respect to the spatial coordinates x a having thereby the following form: 
&(r],x a ) = ^(r),k)e lkaXa and similar for The obtained equations contain 
only gauge invariant variables, so they are valid in an arbitrary gauge. 

Although the metric perturbations are defined in the conformal time 
frame in the sequel the cosmic time t will be used as the function argument 
and all the equations will be formulated with t as the evolution parameter. 

To the background order energy density and pressure are given by fllOH . 
To the perturbed order one has 



-. oo n—1 

o J2 fn J2 (dtSia^dta^-'T + d^Td^- 1 -^) 



n=l 1=0 

ilrrpt nn—X—lr 



n-1 

- 2<5>d t n l Td t n n - L - l T + (19) 
+ 5{U l T)U n - l T + U l T5{U n ~ l T)) - ^{TV(U - Vj nt )5T, 

1 oo n—1 

Sp = - J2 fnJ2 (dtdiU^dtU^-'T + dtU'TdtdiU^-'T) - 

n=l 1=0 

- 2^d t D l Td t n n ' 1 ' l T - (20) 

- 5(U l T)U n - l T - U l T5{U n ' l T)) + ^(TV^ t - V( nt )5T, 

, oo n—1 

tt s = 0. (22) 

where v s gives the perturbed T® components of the stress-energy tensor 
and tt s is the anisotropic stress. Using the Einstein equations one gets that 
tt s = is equivalent to $ = The Bardeen potential \& is proportional to 
the gauge invariant total energy perturbation 

-f + 3(1+ " )ff "1=-i^?*' < 23 > 

The function e is a solution of the following linear differential equation (see 
details in |23j): 

e + H (2 + 3c 2 s -6w)e + 

+ [H(l - 3w) - 15H 2 w + 9H 2 c 2 + )e + -^-A = 0. K ; 
V a 2 J a 2 Q 
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Here w = p/p is the equation of state parameter, (? s = p/p is the speed of 
sound, k = \/k a k a is the comoving wave number and 



a = sp - s 6 + (1 - a V = iL^M £ /n ^ a i n'T5(n"- 1 -'r) - 

k 6 + P „=1 Z=0 

00 n— 1 1 

-E /"E (^^P™"^ + □ / T5(D"-'T)) + -(TC - VU)8T. 

71=1 l=0 9o 

The latter quantity is identically zero for a local scalar field, i.e. in the case 
J-"(D) = f\U + /o- Therefore, A 7^ is the attribute of the non-locality here. 

For the linearized model ([5]) we can consider the background solution 
as given by ([7]) to obtain A in the more convenient form. To do this the 
following relation is useful 

n-1 

<5(D n r) = n n 6r + ^ □ m (5Dp n - 1 -"V. (25) 

m=0 

Using ([7]) and the well-known formula 

n-1 



E •• , r _ 1 

m=0 



one has 

1— in _ jn 

5(a n r) = + £ (SO)*. (26) 

i 

Perturbing the equation of motion for r, one has 

00 

5(Tr)=J2fnS(O n r) = 0. (27) 

n=0 

More explicitly this equation can be written as 

s(tt) =tJ2 (uhr^ 5 ^ + = (28) 

where we have put dr = Yl <5h • 

i 

It follows from ([8]) that if for some Jk we have = as a background 
solution, then <5rfc, contributes only to the second order in the energy- 
momentum tensor perturbations. In this paper we consider perturbations 
only to the first order, and therefore for all Tfc = we can put drt = with- 
out loss of generality. If J- has an infinite number of roots, but we select 
as a background the function r, which includes only a finite number of r^, 
then only a finite number of perturbations dr^ give contribution to the first 
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order perturbation equations, whereas in the second order all perturbations 
are important. 

After some algebra one can get the following expression for A 



A 




m,l 



•F (Jm)J~ l) J mTrnTrnTl Cml j 



(29) 



where £, 




-rf satisfy the following set of equations 




(30) 



Equation (124ft with the above derived A and equations (|30p in a closed form 
describe the perturbations in the case of linearized model. Comparing these 
equation with the equation for perturbations in a system with many fields 
(see eqs. (82) and (85) in [38]) we see they do coincide. Thus perturbations 
become equivalent in the model with one free non-local scalar field and in 
the model with many local scalar fields. In our model, however, the quantity 
which should be considered as energy density perturbation is e. Functions 
Qij play auxiliary role and normally should not be given an interpretation. 

3.2 Curvature and entropy perturbations 

To understand better what is going on in the discussed type of models it 
is instructive to see how the curvature and entropy perturbations behave 
in our model having infinitely many scalar degrees of freedom (interesting 
results in analysing this parameters in case of many fields can be found for 
instance in |39|). 

Comoving curvature perturbations can be expressed as 



where is the Bardeen potential connected in turn with e as (|23|) . Entropy 
perturbations defined as e = 5p — c 2 s 5q can be found as 



Both quantities are gauge invariant and play crucial role in computing vari- 
ous spectral indexes. 

In order to figure out the behavior for the curvature and entropy pertur- 
bations it is enough to find out the behavior of only two functions: e and A. 
Such a generic form in formulae pip and (j32|) persists because the theory 



K = * - -r-(\fr + HV) 



(31) 




(32) 
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can be written in the local form. It is nevertheless obvious that it would 
be difficult to reach similar results in a general case, when Vi n t 7^ and 
equation ([3]) is not linear in T. Moreover, it is not transparent that one can 
make a significant progress for a general operator T even in the linearized 
model. 

There are several specific situations in which we can outline the strategy 
of the further analysis. 

• The first one is already mentioned above and refers to the background 
configuration where only the finite number of scalar fields is excited. 
Then in the linear perturbations we will have no impact of the pertur- 
bation modes related to the zero background fields. This statement is 
obvious having in mind that the local action is quadratic and diagonal 
in scalar fields. 

• The second situation corresponds to particular form of the function T 
such that its roots form a sequence Jj and one can arrange them such 
that Ji — Jj + i —7- when i tends to infinity. Then in the asymptotic 
regime where solutions for the scalar fields are dumped plane waves 
(fl~5l) one gets the factor in the RHS of equation (j30j) tending to zero 
and equations for corresponding Cu+i become homogeneous. In this 
case one can easily solve this equation which now reads 

Cu+i + cCu+i + (-3H + ^ C«+l = 0. (33) 

Here we assumed the coefficient in front of the first derivative to be 
a constant while the background solutions are the planewaves. Two 
modes come out when the time grows 

C«+l ~ Ci + C 2 e- Ct (34) 

and therefore we see that there is a decaying or growing mode depend- 
ing on the sign of c. This constant in the asymptotic regime is given 
by 

c = 3ff -(^+v / %)- (35) 

Recall that the stability of the asymptotic solution requires 3ie( J{) < 
(see (|14p ) and it guarantees that c > and perturbations vanish. Even 
though it is not the final answer and one still has to solve coupled 
equations the claim is that only a finite number of Qj are really coupled 
with e while other functions Qj just produce the inhomogeneity in 
equation for e ()24p . Mathematically we can put all this decoupled 
modes to be zero but than the problem becomes completely equivalent 
to the previous case with a finite number of fields in the play. Here, 
however, we keep trace of other fields even though the are effectively 
decoupled from the system of equations. 
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One example of the operator giving such a behavior is 

F(J) = a-e f3J 

Roots are given by 

J fc = -(log (a) + 2nik) 

Resummation in A (|29|) for all the components which are effectively 
decoupled (assuming this happens from some k = ko) gives 



(36) 



(37) 



Ar 



2Co 



-ct 



where C2 and c are from eq. (|35[) when i = k$. 



(38) 



The third situation is when the values of y/Jl are equidistant meaning 
— V~Ji = cj. Then equation ([30]) in the asymptotic regime 
becomes 

+ (3H + 2^ + cjj Q+ii + ( -•>// + 



> I Ci+li 
7 z ' 



CJ 



E 



Q+P 



Ci+lm Sim ) ~\~ 



(39) 



1 + w 



The advantage with respect to the general situation is that the the 
inhomogeneous part of this equation is universal for any number i 
meaning that we can construct homogeneous equations consdering the 
difference of the latter equations for some i and j. 

It is important to say that all these cases even giving some insight into the 
problem are not very simple to analyze. We hope to address this issue in 
the future analysis of models of this type. 



4 Complex roots J in the linearized model 

4.1 One pair of complex conjugate roots J\ — J and J 2 = J*. 
The background 

If a complex number J is a root of J-, then J* is a root of J- as well. System 
(fTT|) becomes: 



H 2 = -ttGn 
3 



H 



AttG n 
9o 



TV) 



f' 1 ■ .It 1 ) ■ J " iJ '' ' < - ' 



2g a 



T + J*T* ) +A 



(40) 
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In terms of real fields (ft and ip such that t = (ft + iift, t* = <ft — iift,we get the 
following kinetic term: 

£ fc = — ^f 2 + ^^r* 2 = ^U 2 -^ 2 +2^0V, (41) 

2g a *g do v ' 9o 

where d r = 5Re(J r '(J)) and = 9m(J'(J)). In the case di ^ Ej~ has 
a nondiagonal form. To diagonalize kinetic term we make the following 
transformation: 



d r + \ld 2 + d 2 

X = (ft + Cxi) : v= -dcft + ift, Ci = 



In terms of x an d f system (|40p has the following form: 



if 2 = -7rG W 



H = n 2 U - V ) 



(42) 



where J r = Re(J), J m = 9m(J), C = ^j^— So, in the 

case of two complex conjugated roots we get a quintom model (for details of 
quintom models see reviews [4]). 

What is interesting (but not surprising, though) one cannot have non- 
interacting fields passing to the real components. Precisely, fields will be 
quadratically coupled in the Lagrangian. It means that the usual intuition 
about field properties (like signs of coefficients in front the kinetic term or 
the mass term) may not work. 

Following the method outlined in Section 2.3 we find the asymptotic 
solution for the scalar fields with constant H = Hq to be 

r = T+ e a+t + T-e a -\ t* = r; e Q +* + r* e a -* (43) 



where a± = ^Ip- ^-1 i wl — 7^2 ) ■ We assume the first term proportional 

to t + does converge and put r_ = 0. Further we define r + = tq and a + = a. 

The first correction to the constant Hubble parameter and to the scale 
factor in case only decaying modes in r are present gives 



and 



H = H + h = H + h (T 2 + T* 2 ) . (44) 
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4.2 Cosmological perturbations in the neighborhood of the 
solution with complex masses 

Configurations with a single scalar field were widely studied and those ap- 
pearing in the non-local models do not have any distinguished properties. 
Roughly speaking configurations with many scalar fields were explored as 
well but we have here new models featuring complex masses and complex 
coefficients in front of the kinetic terms. As it was stressed above there is no 
problem with this for the physics of our models while properties of such mod- 
els, in particular the cosmological perturbations with such scalar fields were 
not studied in depth. Thus we focus on perturbations in the configuration 
with complex roots J. The simplest case is one pair of complex conjugate 
roots where the background quantities were derived in previous Subsection. 

First we note that the only function Qj is ("12 which we shall denote £. 
Thus, there are only two equations in the system. We focus on the asymp- 
totic regime h <C Hq and after some algebra one arrives to the following 
system of equations 



(46) 



(47) 



(g + p) (c + (3#o + a + a*)C+ (-3H + ^e~ 2Hot ^ c) = 

= (±- J -) ([^V>*V 2 -^)«V]c + 2 5 >), 

e + eH {8 + 3c 2 ) + e (l5H$ + 9i^c 2 + ^e' 2 ^ = 

= 2fc 2 JV)j-'(j> 2 a*V 2 T * 2 / j _ r\ e2[ _ Ho+a+a * )t 

(Q + p)afao^ \a a*) 

where we should use 

H = 2h (r 2 a + r* V 

p + p = F'{J)T 2 a 2 + r{J*)T* 2 a* 2 , 
2 _ P{J)aT 2 (a 2 - J) + F'{J*)a*r* 2 (a* 2 - J*) 

Cs ~ P{J)aT 2 (a 2 + J) + P{J*)a*T* 2 (a* 2 + J*) ' 

The latter system of equations is ready to be solved numerically but in 
order to get some insight in what is going on it is instructive to make some 
assumptions about the value J. This makes some analytic progress possible. 

We recall the SFT origin of the model. Practically this means that values 
of J are determined with the string scales while Hq is expected to be much 
smaller. Therefore, it is natural to assume that \V~J\ ~> Hq. This implies 
a « i\fj . Using the explicit expression for t\ = T$e at rj T$e 1 ^ 1 , representing 
a = x/2+iy/2 and introducing x = ^^rl" 2 !! 7 ! 2 \e xt ( the equations of interest 
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can be written as 



cos(yt)x + 2(v / x 2 + y 2 sin(yt + ab) - xcos(yt))x + 
+ (^cos(yt) (6h Vx 2 + y 2 e x{t - to) sm{y(t - t ) - a b ) + 
—2x\/ x 2 + y 2 sin(yt + ab)) X = —2ye 



cos(yt)e + 2a/ x 2 + y 2 sin(yt — ab)e + 



A- 2 



3a 2 i7n 



+ y 2 cos(yt - a c )e 



2k 2 y 



(48) 



(49) 



l H0 / "o 

where all the constant coefficients are real, e and x are real, to is expected 
to be negative and 



ab = arcsm 



y/x 2 + y 2 ' 



arcsm 



fc 2 



x) 2 + ?/ 2 



The most alarming points of the evolution are yt = ^ + nir where the 
coefficients in front of second derivatives become zero. Numeric integration 
may hit problems at these points if the precision is not very high. In the 
neighborhood of these points one has 



t X - 2 X + 2x X = 2e 
2k 2 

te — 2s — 3Hq£ = — T x- 



(50) 
(51) 



For negative x the solution for e around t = is e = Eq + e\t + . . . meaning 
that these points are not singular for the above system of equations. 

A typical behavior for the function e is dumped oscillations depicted 
in Fig. 1. Such a behavior does not depend on the wavenumber meaning 



o.io- 

0.08- 
0.06- 
0.04- 
0.02- 
epsilon(l) o-f 
-0.02- 
-0.04- 
-0.06- 
-0.08- 



10 12 14 16 18 20 



Figure 1: Typical behavior of the function e. 
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that perturbations with complex conjugate scalar fields do vanish. This is 
different from usual models with real scalar fields where different regimes 
exist and most likely growing modes are present. 

Application of the curvature and entropy perturbation analysis (formulae 
(|3"Tj) and ([32])) shows that both quantities decay exponentially in the presence 
of a pair of complex conjugate roots. Moreover accounting of the effect of 
possible other modes which correspond to other roots with the property that 
these modes are decoupled from the full system of perturbation equations 
(as in (|33]) - (|38|) ) does not change qualitatively the result. More comprehen- 
sive analysis on what happens when complex roots which are significant for 
perturbations are present can be found in [44]. 



5 The linearized model with double roots 

Let us consider an analytical function J- (J), which has simple roots Jj and 
double roots and the function 

Ni N 2 

TB = Y^n + Y^f k , where (□ - Jj n = 0, (□ - J k ff k = 0. (52) 

i=l k=l 

The fourth order differential equation (□ — Jfc)(D — Jk) T k = is equivalent 
to the following system of the second order equations: 

(□ - J k )f k = a k , (□ - J k )a k = 0. 

It is convenient to write \3 m f k in terms of f k and a k : 

U m f k = J k n f k + mJ™- 1 cj k . (53) 

The energy-momentum tensor, which corresponds to tb, has the following 
form |25| : 

Ni N 2 

(r B ) = Y, T ^ (n ) + T ^ fa ) , (54) 

i=l k=l 

where 

T^(t) = \(e^(t) + E^(t) - 9(a> (g p °E prj (T) + W(t)) 
9o 

E,u{n) = ^d,nd„n, w{ n ) = *£S£Lt* , 

Efiuifk) = — ^ (d^fkdyOk + dvfkd^ak) + — d p a k d u a k , 

w( ~, Jj^Vkj - , ( ~JkT"\Jk) , F'(J k ) \ o 
win) = 2 Tk ° k + I — 12 — + — i — J k ' 
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Formulae (|52p and (|54p generalize formulae ([7]) and ([8]) for the functions 
J 7 (J) with both simple, and double roots. Using ([25]) . ([53]) and 



71- 1 

E 

m=0 

we get 



mx 



„m— 1 



'/./• 



n-1 
m=0 



d / 1 — x n \ (n — l)x n — nx n + 1 
" (1 -x) 2 



dx V 1 — x 



n-1 



5{u n T B ) = u n {5r B ) + £ □ m (<yn)n n - 1 - m r B 

m=0 

I — in rn 

iV 2 

+ ^( a "+("- i )f-^r la ) ( ^ ) J ! 



(55) 



(□-J fc ) 2 

To formulate the perturbation equations we note that from (|54p one finds 



N 2 



q{tb) - p(t b ) = 2 E W(n) + X: Wft) 



vi=l 



fc=l 



so that, 



Ni 



5q(t b ) - 5 P (r B ) = 2 £ JiP^ndn + 
i=i 

N'2 r _ 

+ £ J k T" '(J k )(a k Sf k + ffc(5cr fc ) + 

k=l L 



JkF'"(Jk) 
3 



and 



iV 2 

+ E 

fc=i 



a(^+p) 



i=l 



(& k Sf k + f k 5a k ) + 



(56) 



where 



/ ATi iv 2 \ 

e(r B ) +p(r B ) = 2 E £^( ri ) + ^ ^oofo) 

\i=l fc=l / 

Then for one double root we obtain 

A(f fe ) = ^^{^J k P\J k )[a k f k 5a k + alf k Sf kj 
+ 2,F"(J k )a k alSf k + J k T"\J k )o k a k f k ba k 



+ 



(57) 



Note that A(f fe ) / 0, because T"(J k ) / 0. 

After the diagonalizatiorH of the kinetic part of the energy-momentum 
tensor, one can use the general formulae for perturbations in cosmological 



'Explicit formulae are given in [25] 
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models with many scalar fields [38] and get the closed system of equations 
for perturbations. So, we can conclude that in the case of the function -F(D) 
with both simple and double roots we get the system of local equations. We 
plan to consider deeper the case of the function J-(0) with double roots in 
our forthcoming paper. 

6 Summary and outlook 

The main results of this paper are the construction of the perturbation equa- 
tions in non-local models with one scalar field and arbitrary potential and 
consideration of the very intriguing example of perturbations. Namely, using 
the possibility to construct a local equivalent model with many scalar fields 
we find out that masses of these local fields may easily become complex and 
such a case constitutes the above-mentioned example. The characteristic 
feature of the present setup is that all the local fields in fact are not physical 
and play a role of auxiliary functions introduced for the reduction of the 
complicated non-local problem to a known one. As it was noted in [10] lllj. 
for a very wide class of the SFT inspired models the local counterpart is not 
yet studied. Looking strange such configurations do not produce a problem 
for the model since they are not physical quantities. 

Perturbation equations for this local model are (|24p and (|30p where only 
N — 1 functions £ij are independent. The discussion on how a cosmological 
constant can be generated during the tachyon evolution is presented in [71 [TO]. 
We note that perturbations in a quintom model very close to our setup with a 
phantom field without potential and an ordinary scalar field with quadratic 
potential were studied in [42J. Perturbations in models with many scalar 
fields were studied in literature considering various cosmological scenarios 
[38] H3J. 

In the present paper we have worked the indicative example where two 
scalar fields with complex conjugate masses are present. We have demon- 
strated numerically that in the case \VJ\ > H the gauge invariant energy 
density perturbation associated with the matter sector does decay in all 
wavelength regimes in contrary to ordinary scalar field models. The gen- 
eral case of complex masses deserves deeper investigation and is partially 
considered in [44] . 

Moreover we singled out configurations when really an infinite set of 
scalar fields may present but on the other hand there is a chance to analyze in 
full the system of equations for perturbations. These are configurations when 
\/Ji+i — \fJ~i — > and when \J Ji+\ — \fTi = cj = const. This should definitely 
help in understanding how important quantities like curvature perturbations 
and entropy perturbations behave. Also the case of double roots is addressed 
and it is shown that one can again formulate the local equivalent theory and 
build a closed system of equations. However, the deep analysis of all these 
regimes seems to be rather involved and therefore we put it as still open 
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question to be addressed separately. 

Looking further it is interesting to consider perturbations in other non- 
local models coming from the SFT. For instance, models where open and 
closed string modes are non-minimally coupled may be of interest in cosmol- 
ogy. An example of the classical solution is presented in . Furthermore it 
should be possible to extend the formalism presented in this paper to other 
models involving non-localities like modified gravity setups [321 133j . 
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